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Abstract 

Different aspects of the connection between the Bessel process and the conformal 
quantum mechanics (CQM) are discussed. The meaning of the possible generalizations 
of both models is investigated with respect to the other model, including self adjoint 
extension of the CQM. Some other generalizations such as the Bessel process in the wide 
sense and radial Ornstein- Uhlenbeck process are discussed with respect to the underlying 
conformal group structure. 
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1 Introduction 

The Bessel process is one of the building blocks of stochastic processes because of its appli- 
cations and also its simplicity and richness. The Bessel process describes the norm of the 
Brownian motion in arbitrary dimension [1, 2j. This process describes the movement of an 
arbitrary point on the real line in the Schramm-Loewner evolution [3]. It is also important 
in the probabilistic description of some financial markets specially Cox, Ingersoll and Ross 
(CIR) model [4]. Since the Bessel process is connected to the movement of a free random walk 
in arbitrary dimension it could be connected to the quantum mechanics of the free particle. 
The interpretation of quantum mechanics as the classical stochastic equation is a long story 
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and was discussed in [5]. It is based on looking at the generator of the stochastic process 
as the Hamiltonian of the quantum mechanics. Although by now this interpretation is quite 
well known and was discussed also in many books our case, the Bessel process, has not been 
discussed as much as it deserves. The corresponding quantum mechanics is called conformal 
quantum mechanics and it has many applications in different areas of physics such as the 
Calogero model [6j , conformally invariant quantum mechanical models [71 [8] , dynamics of 
quantum particles in the asymptotic near-horizon region of black-holes and dipole-bound 
anions as anisotropic conformal interaction [lOj. Although this quantum mechanics was dis- 
cussed from many different points of view, its connections and similarities with the Bessel 
process have not been addressed so far and that is our goal in this paper. From the stochas- 
tic processes point of view we have a continuous process which is also scale invariant. The 
same interpretation is true in the quantum mechanics part, but the story in the quantum 
mechanics side can be generalized by considering some special boundary conditions at the 
origin which can have stochastic interpretations. It is also possible to define other generalized 
Bessel processes which have interesting well-known quantum mechanical counterparts. We 
will address different aspects of these similarities and connections in this paper. The paper 
is organized as follows: 

In the second section we address many different properties of the Bessel process, the 
starting point is our motivation from the norm of the Brownian motion but we will soon 
generalize the definition in many different directions. In the third section by starting from 
the quantum mechanics of the free particle we will extract conformal quantum mechanics as 
a quantum mechanics with conformal symmetry. Similar to the second section we will soon 
generalize the definition by allowing also a real and even an imaginary number of dimensions. 
Of course in every extension level we should be careful about the physical meanings and 
possible interpretations with respect to the original motivation. In this section we will also 
discuss the self adjoint extension of the corresponding quantum mechanics. In the fourth 
section we will discuss many connections and similarities between Bessel process and CQM, 
in particular we explain the possible meaning of the extension in each side with respect to 
the other. This comparison will also propose the necessity of studying generalized Bessel 
processes. In the fifth section we will investigate the most natural and simple generalizations 
of Bessel process including Bessel process with constant drift, Bessel process in the wide 
sense, CIR model and Morse process which are related to Coulomb potential, free particle 
conformal quantum mechanics, radial harmonic oscillator and Morse potential respectively. 
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We explain how one can get the connections and also see the internal symmetries. Finally in 
the last section we will summarize our results. 

2 Bessel Process 

To define Bessel process we start with the primary motivation for definition of the process. 
It is just the radial part of (5-dimensional Brownian motion, i.e, Rt = \J Bf (t) + ... + i?| (t) . 
Using Ito's formula one can write the following stochastic equation for the Bessel process 

dRt = ^—^dt + dBt (2.1) 

where Bt is the one dimensional Brownian motion. It is easy to see that the above process has 
scaling property, i.e. if Rt is a Bessel process with starting point x then the process c~^Rf.2t 
is also Bessel process with starting point at x/c for positive c. In this level one can consider 
the above equation with the arbitrary real 6. The Bessel process with positive S was studied 
extensively in the literature [H [2] . Before establishing the properties of Bessel process it is 
worth mentioning some properties of Brownian motion in arbitrary dimension. 

It is well known that Brownian motion in (5 < 2 is recurrent and transient for 5 > 2. With 
transient we mean that ii 6 > 2 then almost surly limt_^oo \Bt\ = +oo and the process is 
recurrent if the set {t : Bt €z U} be unbounded for all the sets D in R^; in other words the 
process Xt said to be recurrent if Px{Ty < oo) = 1 for all x,y £ D where Ty is the hitting time, 
the first time at which the process hits point y, and Px{Ty < oo) is the probability of having 
finite hitting time of the point y for the process with starting point at x. Another interesting 
property of Brownian motion is related to the probability of meeting two Brownian paths with 
arbitrary starting points in finite time. Two Brownian paths will meet each other if 5 < 4 
otherwise they will not |TT]. In other words almost surely the probability of intersection of 
two independent Brownian paths in (5 > 4 is zero. 

Let's now summarize some of the properties of the Bessel process. Feller's test, see 
appendix B, indicate that this process has a natural boundary at infinity and a boundary 
at origin which is natural if J > 2, regular if < < 2 and absorbing if 5 < 0. One can 
summarize also the following properties 

/ : for 5 > 2 the Bessel process is transient, 

// : for 5 < 2 the Bessel process is recurrent, 

/// : for 5 > 2 the origin is polar, it will not be touched by the process, 
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IV : for < S < 2, origin could be a killing point or reflecting point. 

The same properties are true for the squared Bessel process defined by Zt = which 
satisfies the following equation 



dZt = 6{dt) + 2^\Zt\dBt. (2.2) 

Using Ito's formula the generator of squared Bessel process is 

Cf{x):=2xf{x)+5fix). (2.3) 

Then one may introduce the Green's function Gx{x, y) as the Laplace transform, with respect 
to time, of the transition density, the density of finding the process starting from x at y, if 
of the process 

Gx{x,y)= e~^'p{t,x,y). (2.4) 
Jo 

The Green's function satisfies the following equation 

CGx{x,y)-XGx{x,y) = 0. (2.5) 
The solution may be factorized as follows 

s \ w^^i^x{x)(l)x{y) iix<y 
Gx{x,y) = < , (2.6) 

[ Wx ip\{y)<P\{x) iix>y 

where wx is the Wronskian wx '■= Tp'^{x)(j)x{x) — ipx{x)(l)'^{x). Since in the case of < 5 < 2 the 
process is touching zero we need to specify the boundary at the origin is a killings boundary 
condition or reflecting boundary condition. Using the above boundary conditions one could 
find the following solutions 

{x^^ Iy{\/2\x) if5>2orO<(5<2 with refiecting origin 
_^ (2-7) 
X 2/_j^(v2Ax) if < or < 5 < 2 with killing origin, 

and 

(t>x{x)=x-'iK^{V2}^), (2.8) 

where 1^ and denote the modified Bessel functions with index v = It is easy to see 
that Wx = \- Using the above solutions one can write the transition density for all 5 > 
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except the killing case as follows 

p{t,x,y) = l(^)-ie-^/,(^) if x>0, (2.9) 

2 _y_ 

p(t,x,y) = -, — , ^, -e 2t if X = 0. (2.10) 

' (2i)i+^r(l + 1/) ^ ^ 

The transition density for < (5 < 2 with the killing origin is 

p{t,x,y) = -(-)-2e-— /_,(V_^). (2.11) 

Similar results could be calculated for the the Bessel process, in this case the generator has 
the following form 

ji:f{x):=^f"{x) + ^{S-l)nx). (2.12) 

Since this process is just the square root of squared Bessel process one can get the transition 
densities in this case by just the transformations x ^ x^ and y ^ y^. For example for the 
reflecting origin we will have 

p(t,x,y) = iy(^)-^e-^/.(f )• (2.13) 

One could get the same answer by the method which is more familiar for physicists and that 
is by using the Fokker-Planck equation which has the following form for the Bessel process 

dtp{t,x,y) = hdl-d,^^^)p{t,x,y). (2.14) 

One can do the transformation p{t, x, y) = x^^'2Q{x, y, t) following by the Laplace transform 
and get the equation 

- IdlQx + '^^^Qx = -AQa- (2.15) 

where Q\ is the Laplace transform of Q{x,y,t). The above equation is just the modified 
Bessel equation with the modified Bessel functions as solutions. It is easy to see that the 
above eigenvalue problem by the change of variable Sx = x'^^'^^/'^^Qx is equivalent to (—5^ — 
(^^^)dx)S{x) = —XS{x) which the operator is the same as the generator of Bessel process 

stated before. 

So far we just addressed the Bessel process with natural boundary conditions in origin 
and infinity but it is also possible to investigate squared Bessel process with (5 > starting 
from the x < or squared Bessel process with S < and arbitrary starting point. These 
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processes were studied in [TJ] and have the fohowing properties: in the case 5 = the process 
will reach zero and stays there. For a squared Bessel process with 5 > and starting point 
X < one could show that it behaves like —ZZ^ process, with starting point at —x, until it 
hits origin and after that it behaves like Zq. Similar relations are valid for the process with 
5 <0 and x > by just reversing the sign of dimension, starting point and the process. For 
example the process with negative dimension and negative starting point behaves as —ZZx- 
It is worth to mention that since for negative dimensions the squared Bessel process become 
negative the square root of it, which is Bessel process, will become purely imaginary and so 
one could be careful about extending the results to the Bessel process, however up to the 
time that the process is positive one could define the Bessel process as well as a real process. 
To complete the discussion we give the transition density of the squared Bessel process with 
negative dimension given in [12] for starting value x > 

P[t,x,y) = k{x,y,6,t)e 2t / ^^exp(— (xu;H )); 

JO w""/^ zt w 

kix,y,S,t) = — ^xi-f|yr(i+i)t^-i. (2.16) 

Another interesting process related to the Bessel process is the time reversed Bessel pro- 
cess. If Rt be a Bessel process with starting point on the positive real line with dimension 
smaller than two then the time reversed Bessel process, defined by (i?(r^^g)„s, s < Tx^o), 
has the same law with a Bessel process Rs,s < Lq^x starting from origin with 6 = A — 5 and 
Lq—,x = {sup t\Rt = x}. In the above notation Tx—ty = \n.i{t\Rt = y} is the first time that 
the process starting at x hits the point y. To make it more clear consider the ensemble of all 
Bessel paths that ended up when they hit the origin in the first time, then consider all these 
paths in the reverse time, the laws of these two processes are the same. The other way to 
say the above statement is as follows: For every bounded function F one can write 

i?o'[i^(i?(T_o)-.,^ < Tx-,o)] = eI-\f{Rs,s< U^x)]. (2.17) 

where F is a function on the set of realizations of the process. The above equality relates 
Bessel process with dimension 5 to Bessel process with dimension 5' = A — 6. One can write 
this duality by defining n = \ ^ 'tj ~ \- Using the equality one could easily see that 

n = sin(^) = sin(^). The above equality is valid for all real values oi 5. It is also possible 
to write the equality as + z) = 0, then one could say that these two processes have the same 
Q{x,y,t) if we work in the positive range of dimensions of Bessel process. In fact we will 
see in section two that these two different Bessel processes are related to the same conformal 
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quantum mechanics. It is also easy to see that 5 = 2 is self dual and 5 = 4 is critical which 
beyond that our reverse time process has negative dimension. The above equality is also 
useful to get the law of certain first hitting times of Bessel process jl2] . 



3 Conformal Quantum Mechanics 



In order to introduce conformal quantum mechanics we follow the same strategy that we 
had in the previous section, we start from the free particle in 6 dimensions. The Schrodinger 
equation for a free particle in radial coordinates has the following form 

1 (3,1) 



2r^~i dr dr 

or after differentiation one could write the following Bessel generator eigenvalue equation 

- ^ - = 2E^ir). (3.2) 

or^ r or 

As it is evident a similar equation was derived for the Bessel process in Laplace space in the 
previous section. By ip{r) = r~^'^^^/^)Q(r) and extending the range of 5 to all real values we 
will have the following Hamiltonian for conformal quantum mechanics 

where p = i-§p and [r,p] = i at the quantum level. The above quantum mechanics as a singular 
quantum mechanics was discussed extensively soon after the discovery of quantum mechanics 
\13\ [T4j and references therein. It is easy also to get the same quantum mechanics for two 
free particles moving in 5 dimensions with ^ interaction by just forgetting the momentum 
of the center of mass term and going to the radial part of the spherical coordinates. The 
conformal quantum mechanics has many interesting properties. Before investigating different 
solutions of the above system we study some symmetries of the system. The Hamiltonian 
at the level of one dimensional quantum mechanics with usual inner product for Q{r) is 
symmetric with respect to — > — but at the level of <f{r) is not symmetric with respect 
to the same transformation. The action corresponding to the Hamiltonian (j3.3p is invariant 
under the following transformations 

^, at + b 



ct + d' 

r'(t') = with ad-bc = 1. (3.4) 

^ ' ct + d ^ ^ 
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The above transformations are the conformal transformations in + 1 dimensions [7]. The 
basic transformations are time translation, dilation and special conformal transformation 
with the following Noether charges 

t' = t + b, H = H = -{p^ + 3.5 

t' = a^t, D = tH - -{rp + pr); (3.6) 

t' = — - — K = t'^H --t{rp + pr) + -r'^. (3.7) 
ct + 1 2 2 

The above generators verify the algebra of the conformal group S0{1, 2) which is 

[H, D] = iH, [D, K] = iK, [H, K] = 2iD. (3.8) 

Using the new definitions Lq = ^{H + K) and L±i = ^{H — K + 2iD) one can write the 
algebra in the more familiar form 

[Lo, L±i] = ^L±i, [i+i, i-i] = 2Lo. (3.9) 

Explicit dependence of K and D on t will guarantee their conservations. The important 
thing to mention is that the argument is not considering the most general case because we 
already knew that the action is invariant up to a total derivative which could be non zero 
in the presence of the boundary in the origin. We will discuss both conformal invariant and 
anomalous case. This will be more clear when we discuss the self adjoint extension of the 
quantum mechanics. Another important thing to mention is although it seems that we found 
three conservation laws for our two-dimensional phase system, it is not difficult to check that 
they are in fact related by the following relation 



HK + KH - 2L>2 = 




(3.10) 



which is also the Casimir operator of the group. In order to present diff'erent aspects of the 
above quantum mechanics we introduce another variable 

9 = ^^-\- (3.11) 

Our quantum mechanics has different properties with respect to the value of g, some of which 
we will summarize in the following. For an arbitrary value of g it is easy to show that if one 
could find one bound state with energy E then scaling the position by an arbitrary factor 
a it is easy to construct a new solution with energy a'^E, which means that if there exist 
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any bound states then there is a bound state for every negative energy. This is a direct 
consequence of the scahng property of this model. The same story is true for the positive 
energy solutions of the model which means that we have just planar waves with all the possible 
positive energies which was also used to obtain equation (j2.15p . In the previous section we 
considered lim,._^o '/'(^) = for the killing boundary condition and for the reflecting case 
we had lim^^o ^^[r^ =0. In the level of conformally invariant quantum mechanics we will 
ignore the negative energy solutions which means that our Hilbert space is made by wave 
functions corresponding to continuous positive energy solutions. The corresponding boundary 
condition also will be discussed in the end of this section. 

The Hamiltonian could have other solutions dependent on the value of the coupling g. 
Firstly we should mention that the Laplacian operator that was the starting point is self 
adjoint if we consider it in the whole space. But one could extract other solutions by removing 
the origin from the domain and considering the self adjoint extension of the operator in the 
new domain. Of course existence of the extension is dependent on the value of g or in other 
words to the corresponding dimension. All the necessary aspects of theory of self adjoint 
extension were discussed in the appendix A. In the case of our quantum mechanics the 
possible extensions were discussed in [T5l [T6| [T7]. 

For g > J the Hamiltonian is self adjoint with removed origin now and has a scattering 
sector with the following solutions 

ip{r) = {V2Er)^ Ju{V2Er) or {V2Er)^Y^{V2Er), (3.12) 

where J and Y are the Bessel functions. For > 1 just J can be considered because in 
this case Y at the origin is infinity and for < — 1 just Y has desired property. For i' > 1 
which is equivalent to 5 > 4 one could argue that it is not possible to have a bound state 
solution because always we can just consider one of the solutions of the Bessel equation. The 
boundary condition in this case is lim,._+o = and can not be extended. 

For g < J since the Hamiltonian is not self adjoint one can find a required self adjoint 
extension. It is better to distinguish between the domain ^ < 5 < | and g < First we 
discuss the domain < g < \ which is equivalent to — 1 < < 1 or < < 4. In this range 
the Hamiltonian requires a self adjoint extension with the self adjoint parameter z which is 
responsible to map two deficiency subspaces by the unitary map e*^. The important point to 
mention is that although the Hamiltonian in the domain D{H) = {95(0) = ^''(0) = 0} is not 
self adjoint, it is still Hermitian. 
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If we consider as the adjoint of H, with the same differential representation as H, then 
from the Von Neumann's theory of deficiency indices we know that the deficiency subspaces 
K± are made by the square integrable solutions of the equation H*(j)± = ^i(t)± in the desired 
domain. In our case the solutions are 

(/>+(r) = AHl{re'^), (3.13) 
0_(r) = AHl{re-'^), (3.14) 

where H}^ is the Henkel function and both of 4>±{r) are square integrable in the half line. 
Using the above solutions one could argue that the Hamiltonian is self adjoint in the domain 
Dz{H) = D{H) © {n((/)+(r) + e*^0_(r))} where u is an arbitrary complex number. To 
find the valid boundary condition we need ifj to be in the self adjoint domain. Consider 
$ = <l)+{r) + e*^0_(r) then ip is in the self adjoint domain if < ^\Hip > = < H^lip > or 

Iim[cl>*^-V;^] =0. (3.15) 

r^o^ dr dr ' ^ ^ 

To check this we need the behavior of $ close to the origin 



where 



^(r) 7 — r[Ar''+2 + Br~''+2]^ (3.16) 

smfvrz^ 



2'^r(l + I/) 



B = , , , (3.17) 

2-^r(i-i/) ' ^ ^ 

for 7^ 0. Using the above relations one can write the equation (j3.16p for the boundary 
condition as follows 

(^^-+1/2 + ij^--+i/2)^ _ (^(^ + ^)r-"-^/' + Bi^- - v)r-^-^'^)xl, ^ 0. (3.18) 

The equation for = ^ is like B^\q = A%l)(fi). Actually = | is very interesting because 
firstly it corresponds to 5 = 3 and also because it describes the possible boundary conditions 
for the free particle in the half line. The simplicity of the results helps us to investigate the 
possible meanings of the above self adjoint extension with respect to the Brownian motion 
and the Bessel process which we will discuss in more detail in the next section. The next 
important case which deserves separate calculations is = 0. In this case the function <I> has 
different properties near the origin. Using (/>+(r) = 01 (r) and 

2i 1 1 2i 1 • 

<I>(r) — r2 ln(r) + [- H (7 - In2)]r2 + e*^cc 

vr 2 vr 

= {A + A*e'^)A ln(r) + {B + B*e'^)r^. (3.19) 
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where 7 is the Euler constant and cc is the conjugate of the first term. One can find the 
following boundary condition 

Ar-2 Hr)-^ - ^[A(l + — + 1) + ^ 0, (3.20) 

where 

A = A + A*e" B = B + B*e'\ (3.21) 

One can extend the results also for the case with g < —\ which is related to pure complex v. 
In this case we have still a one parameter self adjoint extension and all of the results are still 
valid if we put imaginary in to the formulas, in particular in to the (j3.18p . The important 
point to mention is in fact that by introducing the self adjoint parameter we are classifying 
all the possible physical boundary conditions of our system. In our case we should mention 
that after the self adjoint extension our Hamiltonian is not necessarily scale invariant and 
so could have a negative ground state which is in fact the case here. Of course the ground 
state depends on the self adjoint extension parameter. The energy and wave function for 
our Hamiltonian were discussed before in [16j and references therein, and the results are as 
follows. 

The Hamiltonian does not have any bound state for 5 > | but admits one bound state 
in the range ^ < <? < |. The wave function of the bound state up to the normalization 
constant has the following form 

V'(r) ~ AK^{^-2E^r), (3.22) 

with energies 

sun^2 ' 4 ) 

E,=o = -exp|cot(^). (3.24) 

For the case g < ^ we have infinite discrete bound states |13j with the same wavefunction 
as above but with imaginary z^. To get the energies we just need to put the wave function in 
the corresponding boundary condition which yields 

A 

E = -2{—yexp{2aTgT{l + u)-2mT}, for n G Z. (3.25) 
B 

It is worth to emphasize again that since after the self adjoint extension the action is not scale 
invariant we do not need to worry about having a bound state because it does not enforce 
any other bound state by the scaling argument. 
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One could summarize this section as follows: we introduced conformal quantum mechanics 
as the quantum mechanics of the free particle in 6 dimensions. Constraining the domain of 
the quantum mechanics to the space without origin forces new boundary conditions to the 
system. By these boundary conditions for some values of g we need to find a self adjoint 
extension of the Hamiltonian. The method of extension helps to find all the relevant boundary 
conditions and so the physics of the model is related drastically to the self adjoint extension 
parameter. The extended Hamiltonian which is not scale invariant admits discrete bound 
states. Since after the self adjoint extension our conformal symmetry has some anomalies, 
may be the expression conformal quantum mechanics is a bit confusing but because of our first 
motivation we will keep using it. Different aspects of this anomaly were discussed extensively 
in |18l [TUl [TUl fIU\ from the re-normalization group point of view. The above results could 
be also explained from the framework of two free particles moving in 5 dimensions. The 
self adjoint extension is just giving boundary conditions corresponding to the time that two 
particles meet each other. For example, bound states of the above system could be seen as 
the bound states of two particles which behave like one particle. In the next section we give 
some hints on how one can translate the above results in the Brownian motion language. 

4 Bessel Process and CQM 

In this section we are interested in discussing different features of the similarity between 
Brownian motion in 5 dimension and free particles from a quantum mechanics point of view 
and also similarities between Bessel process and CQM. 

The first remark is for the Bessel process with positive dimension. We had origin and 
infinity as the natural boundary conditions which is also the case for our CQM. In addition 
before extending the Hamiltonian we have a scale invariant system which is the case also in 
the Bessel process. Just as we remarked in the previous section, keeping the scaling symmetry 
of the CQM gives us an unbounded continuous spectrum which is also true if we try to find 
the transition density of the Bessel process. To make more clear the connection between two 
models one can use the following relation between the transition density of the Bessel process 
and the path integral of CQM 



r<iii)=y 



P{x,y,t) = / [dq{T)]exp[-S{q{T))] 



Jq{to)=x 
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S{q) = Jl(<i+^-^)'dt (4.2) 

It is not difficult to see that by a transformation p{t,x,y) = x'^^^Q{x,y,t) one can go from 
tlie above action to tlie familiar conformal action that we discussed in the previous section. 

By the above discussion it seems that the self adjoint extension is something beyond our 
normal understanding of the Bessel process. It is not difficult to see that it is related to 
the different possible boundary conditions for the Bessel process at the origin. Different self 
adjoint extensions are related to different boundary conditions or to different measuring of the 
paths of the Bessel process. To make it more rigorous, by Bessel process we mean some kind 
of generalized Bessel process which is related to the diffusion of Brownian motion in the space 
without the origin. From the two Brownian motions point of view naively one could argue 
that since we have a self adjoint extension just for 5 < 4 it is just a hint that two Brownian 
paths could meet in just 5 < 4 dimensions. Of course we are aware that these two statements 
are not equal, but they are related because in the case of two Brownian paths one could 
think about moving one Brownian motion in the space which remains after removing from 
the original space the region filled by another Brownian motion and the boundary condition 
on the fixed Brownian motion is the same as the one that we discussed before. In fact the 
interaction of the paths is a point like interaction and as we will explain soon this interaction 
plays no role in dimension greater than three in quantum mechanics, neither for attractive 
nor for repulsive interactions. 

To understand what could be the meaning of self adjoint extension in the case of Brownian 
motion, it is worth to understand the situation first in the case of quantum mechanics. One 
way to get a motivation for the self adjoint extension in the pointed space is by considering 
the delta function potential at the origin. Of course this will work just for the integer 6. 
The cases 5 = 1,2,3 are studied with more details for two reasons, firstly they have more 
applications and secondly the 5 > 4 cases are trivial in some senses that we will discuss 
soon. The equality of a Hamiltonian with a delta function potential with a free Hamiltonian 
on a space with one point deleted plus a boundary condition were discussed by Jakiew |21j . 
The equality is in the level of re-normalized delta function which gives the same scattering 
data and bound states as produced by the self adjoint extension. One could re-normalize the 
delta function with different functions including sphere delta function, square well and lattice 
regularization. For example in the sphere delta function case, to regularize V{r) = v6{r) we 
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can use 

V{r) = ,^^^5{r-R)Mt\^ C2 = 1, c^ = \; (4.3) 

for dimensions d = 2,3. The corresponding limit is i? ^ 0. In the case of the square well 
we have V{r) = —DQ{R — r), where G is the Heaviside function, D is a coonstant and 
the limit is the same. Since these regularization procedures were discussed extensively, we 
ignore the details but we explain the strategy, for more details see |21| . Using the above 
regularized delta potentials firstly we should solve the problem for the regularized potential 
by considering the continuous wave function at r = ii. This matching leads to equations 
between different parameters including v, R and energy. Then by defining the re-normalized 
coupling V as a function of R and v, one could derive the energy and wave function of the 
system exactly as we found in the previous section. For example for the three dimensional 
case we have by definition i = i + and after matching the wave functions we will have the 
following equations for the bound state energy Ej, and the phase shift of the s-wave scattering 
sector 

V2m> = ^, 4.4 

V 

tan((5o) = -J. (4.5) 
zvr 

Similar results for the two dimensional delta function were given in [21j . 

This regularization procedure will not work for dimensions d > 4 because it is not possible 
to absorb all of the divergent terms in this case. Although it seems that the re-normalized 
attractive delta function in two dimensions has the same physics as the free particle quantum 
mechanics in the plane without origin, this matching is not complete for some reasons. For 
explanation we need to make some remarks on the general properties of the delta function 
potential in arbitrary dimensions. Firstly the repulsive delta function is trivial in dimensions 
higher than one which means that the phase shift scattering is zero for this case and so the 
scattering matrix is equal to one P2j- For the attractive case as we remarked before it is 
not possible to define a zero range potential in more than three dimensions possessing bound 
states with finite energy. For simplicity we focus on the three dimensional case and compare 
with more detail the delta function potential and the self adjoint extension counterpart. Let's 
take the subset of the possible boundary conditions with the following property for each real 
11, 

m^-^t^„. (4.) 
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For positive energy E = the solution is as follows 



V'(^) = -(sin(/cr) + tan((5o) cos(/cr)) 



(4.7) 



where 5q is the phase shift corresponding to the s-wave. It is not difficult to see that for the 
attractive delta function using equation (j4.6p and (j4.7p one could get rj = v. This matching 
can not be done for the repulsive case because we already know that the repulsive delta 
function is trivial, but the self adjoint extension with r] < still has a phase shift and so 
it is not trivial. At the level of the attractive case one may still interpret the physics of 
the self adjoint extension for the integer 6 as the stochastic process in the presence of the 
regularized delta function. However , writing a stochastic equation for the radial Brownian 
motion with one removed point is not obvious and needs more investigation. Of course the 
definition of a stochastic process for the generic case with arbitrary 5 or is more difficult. 
Although the boundary condition for 6 = 3 oi u = ^, which is equivalent to the half line free 
quantum particle, is very simple, it is enlightening. One could see the boundary condition 
()4.6p from the stochastic process point of view at one dimension as follows: r/ = is the 
Dirichlet boundary condition which is an absorbing boundary condition. The particle will 
be absorbed by the origin after hitting that. This also corresponds to the conventional free 
Hamiltonian with scale invariance. From the two particle point of view it is like absorption 
of one particle by the other one when they touch each other. The Dirichlet boundary case 
is just reminiscent of the equality of the 3 dimensional Bessel process with the Brownian 
motion on the half line with absorbing boundary condition. 

The other extreme case is the Neumann boundary condition r/ — > — oo which is a reflecting 
boundary condition at the origin for the Brownian motion. Other negative values of rj corre- 
spond to a mixing of Dirichlet and Neumann boundary conditions. This correspondence was 
discussed in detail in [24:\ [25t [26] and the Green function has the following form for arbitrary 
value of the self adjoint extension 



Gn{x,y,t) = GF{x-y,t) + GF{x + y,t) + — dwe ^""GFix + y + w) r/ < 0; (4.8) 



Gn{x,y,t) 



+ 




(4.9) 




V 

_J_J{x-y)^/2t 



(4.10) 
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For the special cases, Dirichlet and Neumann the results ctrG clS follows 



Grj=o{x, y, t) = Gf{x -y,t) + Gf{x + y, t); 



(4.11) 



(x, y, t) = Gf{x -y,t) - Gf{x + y, t). 



(4.12) 



One can use the above equations to get the Green's function of the Brownian motion in three 
dimensions when the origin is removed. Of course one could get the Green's function for the 
general case by using an orthogonal eigenvalue expansion for the self adjoint operator. The 
Green's function with respect to the solutions of the Hamiltonian has the following form 



In the case oi g < j since we have also a bound state we need to add new terms coming from 
the discrete contribution of this states by adding X]fe(^*'^''*V'(^)V'*(^)) to the Green's function, 
where Ef^ is the bound state energy. The important conclusion is that one could also consider 
the above Green's function as sort of a generalized Bessel process in 5 dimensions. In fact 
one can get the transition density for the above processes by just Wick rotation. The above 
generalization of Bessel process has not been appeared in the mathematical literature. From 
the stochastic process point of view one can derive the above solutions by considering the 
local time of the process [I]. The definition of the local time of the path lo at the point a is 
as follows 



where lx-\-e{Bs) is the indicator for the time that the process is in the interval [0, x + e]. 
One could naively write the above equation as an integral over a delta function as ti{uj,0) = 
Jq 6{x{t))dt. Using the above equation the extended transition density for 5 = 3 is just the 
expectation value of exp (^^|^tj). The corresponding stochastic process for the free particle in 
the half line is called elastic Brownian motion [Ij . To the best of our knowledge the problem 
has not been discussed by the mathematicians for the generic case. For the Bessel process 
that we discussed in the second section one just needs to consider i? = in the equation 
(j3.18p . In this case we will recover conformal symmetry for our process again. It will be 
also interesting if one could get the results for the Bessel process with negative dimension by 
using the quantum mechanics of a free particle, in particular the equation ()2.16p . 




(4.13) 




(4.14) 
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5 Generalization of Bessel Process and CQM 



In this section we want to present some possible generalizations of the Bessel process with 
the well known quantum mechanical counterparts such as, Bessel process with constant drift, 
Bessel process in the wide sense, Cox, Ingersoll and Ross (CIR) model and Morse process 
which are related to a Coulomb potential, conformal quantum mechanics, radial harmonic 
oscillator and Morse potential respectively. Of course not all of the above processes admit 
conformal symmetry; they are reasonable perturbations of the conformal invariant case. 

5.1 Bessel process with constant drift and Coulomb potential 

The definition of the Bessel process with constant drift is given by the following stochastic 
process 

dR^ = {^--^ + fi)dt + dBt. (5.1) 
ZKt 

From the Feller classification the boundary at the origin has the same classification as for the 
standard Bessel process, the infinity is a natural boundary condition which is attracting if fj, 
is positive and non-attracting for a negative drift. The generator of the process is 

Cf{x) := ^fix) + {^{6-l)+ ^,)f{x). (5.2) 

One could map the solution of the equation (j2.5p to the Schrodinger equation of the Coulomb 
potential. To do so we need the Liouville transformation 

1 r s-i EL. r 
i?,(r) = 24(— ) — ev^G,(^), (5.3) 

where r = \/2x. Using the above transformation one could write the equation (|5.2p as 

d'^Rx ,c /(/ + 1),„ ^„ 

^ + (- - ^7^)^A = -ERx. (5.4) 

where E = —A + / = and c = —fJ-^^- For integer 6 the equation is related to the 
energy levels of a Hydrogen atom. The different aspects of the equation (j5.4p were discussed 
in [27| including the transition density for the different values of fi and 6, here we just mention 
some properties. The Bessel process with constant drift admits a discrete spectrum for all 
negative values of /i which is expectable for those people who are familiar with the Hydrogen 
atom. For < 5 < 1 the process admits just one trivial bound state but for 5 > 1 it has 
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infinitely many bound states. The process also admits discrete spectrum for the positive 
values of fx for dimensions in the range < 6 < 1. In this case = 1 is the critical dimension 
for admitting discrete energies or not. 

5.2 Bessel Process in the Wide sense 

Another very important generalization of the Bessel process is by considering a Brownian 
motion in 5 dimensions with a drift of magnitude /x > and then taking the radial part as a 
time homogeneous diffusion in [0, oo). It was shown in [29j that the generator of this process 
is 

cm ■■= Ifir) + - 1) + h-\f,r)^^^)nr), (5.5) 

where 

K{x) = i-mi + u)Ux), iy = l-i. (5.6) 

The transition density of this process is related to the transition density of the Bessel process 
ps(t,x,y) as follows 

PS,fM(t,x,y) = e~''^^^^h;;^{fix)ps{t,x,y)K{ny). (5.7) 

The above process has an interesting time inversion property. It was shown in [28j, see also 
|30j that if Xt with t > is a diffusion process and tXi with t > is homogeneous and 
conservative, no killing in the interior domain, then both processes are necessarily Bessel 
processes in the wide sense up to a time scaled re-parametrization. It is easy to generalize 
the above statement to the processes with time inversion property of degree a by demanding 

a Markov process with homogenous t"Xi. They can be obtained by just considering Bessel 

t 

processes in the wide sense with appropriate power. The Hamiltonian corresponding to this 
process after transformation ps^^{t, x,y) = r'^^^^'^hp{iir)Qp{iir) is as follows 

H = \{P^ + V[r)) 

V{r) = "^^^ + ^1^ (5.8) 

where we used the equality fj? = + A {^h+u{Mr)+^^j^l2+i,{i-ir)) ^ gjj^^g ^j^g third term does not 
have any singularity it means that the general aspects of the above quantum mechanics are 
the same as the Bessel process, in particular it will have a self adjoint extension for the range 
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V <\. This equality means that basicahy the above two distributions have similar spectrum 
and also it is a very simple proof for the equation (j5.7p . One could use all the previous 
solutions for the self adjoint Hamiltonian to get the solutions in this case. The consistency of 
the inner product is coming from the Doob's /i-transform. For the case of Brownian motion 
with negative drift one just needs to replace the modified Bessel function Ii,{p,r) with 
K^{pr). The same as the Bessel process case here also one can generalize the Bessel process 
in the wide sense by considering all the possible self adjoint solutions of the Hamiltonian [?] 
and then using the Doob's /i-transform. It will be really interesting to study this generalized 
cases with respect to the time inversion. 

5.3 CIR Model and Radial Harmonic Oscillator 

The Cox, IngersoU and Ross (CIR) model [1] or radial Ornstein- Uhlenbeck process is widely 
used for interest rate framework such as some stochastic volatility models [31] • The definition 
of the CIR family of diffusions is by the following equation 

dNt = (a - bNt)dt + c.fiN^\dBt, (5.9) 

with A'^o ^ 0, a > 0, c > and b is an arbitrary real number. It is not difficult to see that 
a = 5, b = and c = 2 is the squared Bessel process. A CIR process can be represented in 
terms of the squared Bessel process as follows 

Nt = e-''Z{^^{e'' - 1)), (5.10) 

where Zt denotes the squared Bessel process of dimension 5 = The above equation can 
be checked by using Ito's formula for deterministic time transformation of a squared Bessel 
process. Using the above connection one can easily classify the properties of the origin as a 
boundary for the process; we just need to consider |^ as the dimension of the corresponding 
Bessel process. To make the connection with the radial harmonic oscillator let's first map 
the CIR process to the square root of it Mt = \/\Nt\ by using the Ito's formula as follows 

dMt = ((^ - j)^ - \M)dt + ^dBt. (5.11) 

One can easily map the above process to the radial harmonic oscillator by the following 
Hamiltonian 

H = l\p' + '^x' + ^ + e]. (5.12) 
2 4 
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where uP' = b"^, k = ^(a — — ^ — 1) and e = — |(o — ^) — b. The above Hamiltonian is 
just the radial part of the harmonic oscillator in 6 dimension. The above calculation shows 
that by deterministic time change one can go from scale invariant Bessel process to non- 
scale invariant process with stationary solution. This is a hint to believe that it should be 
possible to move from a conformal quantum mechanics to a radial harmonic oscillator by 
time translation. This is in fact natural and was done long time ago in [7J. The strategy is 
as follows: Firstly one could write S0{2, 1) generators in a more familiar form by defining 
new generators as follows 

5 = b^K-aH), 
2 a 

R = h-K + aH), (5.13) 
2 a 

where a is a constant and the commutators are as follows 

[D, R] = iS, [S, R] = -iD, [S, D] = -iD. (5.14) 

The important point is that the operators D and S correspond to hyperbolic non-compact 
transformations and R is the generator corresponding to a compact rotation. Since all of the 
above operators are the invariants of the action one could define a generic operator as 

G = uH + vD + wK, (5.15) 

as a constant of the motion. Of course it will correspond to compact rotation in three 
dimensions if we consider A = v'^ — Auw < 0. This will be important to get a theory with 
reasonable time evolution. From now on we will just consider this case. The simplest example 
of this kind of operators is R with A = — 1. The action of the operator G on a wave function 
is as follows 

G I ^{t) >= i{u + vt + wt'^)-^ I ^(i) > (5.16) 

which by time transformation could be written as 

G I ^(t) > = ^(r) > (5.17) 
ar 

Aw , ,2wt + v. / V -loN 

7 = —;= ] arctani — — arctan — —)\. (5.18) 

^/Za ^ -A ^ ^-A^^ ^ ^ 

In the new parametrization one could think about G as the new time translation and thus 

the new Hamiltonian. After some algebra one could write the following Hamiltonian 

H = l[p'-j-' + j,l (5.19) 
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as the most general possible Hamiltonian that could be extracted from the 50(2, 1) group 
with applicable time translation. To compare with the CIR model one could write A = —6^ 
and g = k which is an indication to believe that the corresponding time translation in a Bessel 
process is related to a compact rotation in a three dimensional space with metric (—1, —1, 1). 
The generators of rotations in the planes xy, yz and zx are R, D and S respectively. Different 
aspects of the self adjoint extension of the above Hamiltonian were discussed in [32j and 
references therein and they are quite similar to the case without harmonic potential. 

5.4 Morse Potential 

In this subsection for the sake of completeness we want to discuss briefly another related 
physical model, Schrodinger equation with Morse potential [33]. This potential is exactly 
solvable and has many applications in molecular physics [331 [M] . This system is related to 
the quantum mechanics of the radial Harmonic oscillator. One can find the Morse potential 
by just making the variable change u = — 21nx in the Schrodinger equation, Hip{x) = E'ip{x) 
with Hamiltonian (j5.12p . Then we will have the following Schrodinger equation with the 
corresponding Hamiltonian H"^(f){u) = E"^(f){u) where 

ffrn ^ 1 2 ^g-2n _ - e/2 » = --{- + -). (5.20) 

2^ 32 4 ' 4^8 2^ ^ ^ 

The above equality means that the Morse potential is just the canonical transformation, 
i.e. u = — 21nx and pu = —^xpx, of the radial Harmonic oscillator and so the dynamical 
symmetry group of the system is still 50(2, 1) with the transformed coordinates [35] . The 
other important issue is that in this case the domain of the quantum particle is all of the real 
line and so it is not necessary to worry about the boundary condition at the origin. It is also 
easy to extract the corresponding stochastic process by using the Hamiltonian as follows 

dUt = ife-^' + l)dt + dBt, (5.21) 

where f = —j and I = -^{E — e/2) + ^. It is also possible to extract the above equation by 
using equation (j5.1ip and the Ito's formula plus time re-parametrization. The connection of 
the Morse potential to functionals of the Brownian motion was discussed before in [36j. It is 
easy to see by the Feynman-Kac formula that the Kernel of the Morse Potential is equal to 
the following expectation in a stochastic process 

E[expiXkat-h^At)\Bt = y]; (5.22) 
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with A = ^ and k = ^ ^^'^ and 

at= f exp{Bs)ds; At = f exp{2Bs)ds. (5.23) 

JO Jo 

In [36] the connection of the above process to the Maass Laplacian [37] were also discussed 
extensively. 

6 Conclusion 

In this paper we explained many aspects of the connection between the Bessel process and 
its possible generalizations on the one side, and the conformal quantum mechanics and its 
generalizations on the other side. The Bessel process as the path integral interpretation 
of conformal quantum mechanics has conformal symmetry before considering the non-Feller 
boundary conditions of the process which correspond to the self adjoint extension of the 
corresponding quantum mechanics. These boundary conditions could be Feller type if we 
consider more generalized stochastic equations. We also discussed some generalizations of 
the Bessel process that have interesting well-known quantum system counterparts. These 
generalizations are based on the connection between the Green's function of the quantum 
particle and transition density of the Bessel process. Of course there are also many other 
systems but we focused on those that have conformal symmetry as the dynamical symmetry 
of system. This work could be extended in many directions including a rigorous study of the 
stochastic processes that correspond to the self adjoint extension of the singular quantum 
mechanics in the finite domain. This could be done by using the definition of local time. 
Investigating time inversion properties of the generalized Bessel processes in the wide sense 
can be useful in classification of time invertible processes. 

Another interesting study could be the study of the process and quantum mechanics as 
a system with supersymmetry; in this case we will have super conformal symmetry as the 
symmetry of the quantum mechanics. 
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7 Appendix A: Self Adjoint Extension of the Hamiltonian 



In this appendix we will summarize Von Neumann- Weyl method of self adjoint extension 
for the Hamiltonian operators \15\ [38] . 

Consider a Hilbert space TC then an operator {A, 'D{A)) defined on Ti is said to be densely 
defined if the subset T^iA) is dense in Ti., i.e., that for any ip £ TC one can find in 'D(A) a 
sequence (pn which converges in norm to tp, in other words we should have /q°° \ ^/J — (pn]"^ dx < e 
for arbitrary positive e. 

The adjoint operator of an operator H with dense domain 'D{H'^) is H'^ . The domain 
'D{H'^) is the space of functions ip such that the linear form — > {il),H(/)) is continuous for 
the norm of Ti which guaranties the existence of a t(^^ Ti such that 

{i;,H(t>) = {^l^l(P). (7.1) 

Then one may define H^ip = ip^ . An operator {H,T?{H)) is said to be symmetric or 
Hermitian if for (j)^ & ^(H) we have {(j)^Htp) = [Hcp^ip). The operator H with the dense 
domain 'D{H) is said to be self-adjoint if ^(i^t) = 'D{H) and W = H. 

Definition of the deficiency subspaces K,± are by 

/C± = {V- G V{H^), = ±ii)}, (7.2) 

with dimensions n± which are called the deficiency indices of the operator H and will be 
denoted by the ordered pair (n-|_,n_). The following theorem, discovered by Weyl and gen- 
eralized by Von Neumann, is the most important result of this appendix. 

Theorem: For an operator H with deficiency indices (n+, n_) there are three possibilities: 
1: If n+ = n_ = 0, then H is self-adjoint. 

2: If n+ = = n, then H has infinitely many selfadjoint extensions, parameterized by 
a unitary n x n matrix with real parameters. 

3: If n+ 7^ n_ , then H has no self-adjoint extension. 

A relevant example for the above theorem is the Bessel operator discussed in the paper. 
The case on the half line was discussed extensively in section 2 and we will not discuss 
it again but the case in the the finite interval [0, L] is more complicated and needs to be 
discussed separately. In this case for v > 1 only one solution is possible and the deficiency 
indices are (1,1), while for < < 1 both solutions are acceptable and the deficiency 
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indices are (2,2). This is quite natural because for the Bessel operator in the finite interval 
we have another boundary which is like infinite well, so we can expect that the boundary 
condition is like the boundary condition of S = 3 case for u > 1. For the case < < 1 
the situation is more subtle because it is also possible to have some interacting boundary 
conditions. In almost all applications the conditions are separated and so one can look at 
the non-singular boundary, i.e. L, as an infinite well with one parameter extension as the 
5 = 3 case. The boundary condition at the origin is just as before. The interesting point is, 
since for a differential operator of order n with deficiency indices (n, n), all of its self-adjoint 
extensions have discrete spectrums one could argue that for this case all of the energy levels 
are discrete. In the other words for the particle in the finite sphere with origin removed the 
spectrum of energy is completely discrete for < 5 < 4. 

8 Appendix B: Boundary Conditions for the Stochastic Equa- 
tions 



In this appendix we would like to summarize Feller's classification of possible bound- 
ary conditions for the one dimensional stochastic equation |39^ HQ] . Consider the following 
equation as our stochastic equation 

dxt = fi{xt)dt + a{xt)dBt (8-1) 

To classify the possible boundary conditions we need to define the following two functions as 
the scale function s{x) and speed measure m{x) as follows 

Using the above functions one can define the following four different functions for diffusion 
in the interval with endpoints / and r 
ry 

S[x,y\ = I s{z)dz, S{l,y] = lim S[x,y], S[x,r) = lim S[x,y], (8.3) 
M{c,d) = / m{x)dx, M {I, y] = lim M[x,y], M[x,r) = lim M[x,y], (8.4) 

Jc x^l+ x^l~ 

^(0 = r S{l,y]m{z)dz, = r S[z,r)m{z)dz, (8.5) 

Jl Jx 

N{1) = S[z,x\m{z)dz, N{r)= S[x,z\m{z)dz. (8.6) 
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The boundary classification depends on the behavior of the above functions and one can put 
the possible boundary conditions in one of the following four types for the endpoint e: 

1: regular if X^(e) and iV(e) be finite, 

2: exit if Z^(e) be finite and N{e) be infinity, 

3: entrance if J2{^) be infinite and the N be finite, 

4: natural if J2{^) ^^"^ ^i^) t>c infinity. 

For entrance, exit and natural, no boundary conditions are needed but for regular bound- 
ary, the conditional probability is not unique and dependent on the boundary conditions. An 
exit boundary can be reached from the interior point of the domain with positive probability 
however it is not possible to start the process from the exit boundary. An entrance boundary 
cannot be reached from the interior point of the domain but it is possible to start the pro- 
cess from the entrance boundary. A natural boundary cannot be reached in finite time from 
the interior point of the domain and it is impossible to start the process from the natural 
boundary. A regular boundary is accessible and could be reflecting if m{e) = and sticky 
if m(e) > 0. For reflecting boundary point the process spends no time in that point but for 
sticky boundary point the process spends a positive amount of time at a sticky point. For 
example m(e) = oo is called a killing boundary condition. 

For the Bessel process the scale function and speed measure are given by 



Different boundary possibilities for the Bessel process were discussed already in the paper. 
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